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We describe an efficient DQCl-algorittim to quantify tlie amount of Geometric Quantum Discord 
present in the output state of a DQCl computation. DQCl is a model of computation that utilizes 
separable states to solve a problem with no known efficient classical algorithm and is known to 
contain quantum correlations as measured by the discord. For the general case of a (1 + n)-qubit 
DQCl-state we provide an analytical expression for the Geometric Quantum Discord and find that 
its typical (and maximum) value decreases exponentially with n. This is in contrast to the standard 
Quantum Discord whose value for typical DQCl-states is known to be independent of n. We 
experimentally demonstrate the proposed algorithm on a four-qubit liquid-state nuclear magnetic 
resonance quantum information processor. In the special case of a two-qubit DQCl model, we 
also provide an expression for the Quantum Discord that only requires the outcome of the DQCl 
algorithm. 

PACS numbers: 03.67.-a, 03.67.Lx, 76.60.-k. 



Since its inception in 2001, Quantum Discord (QD) [1, 
[2] and related measures have been used to quantify the 
amount of non-classical correlations in a physical system. 
In addition to quantifying the most well known quantum 
correlations of entanglement, they also capture the quan- 
tum correlations that exist in separable states. Discord 
measures are based on the premise that if a measurement 
on one part of a bipartite state disturbs the total state, 
there must be correlations stronger than what is found 
in the classical world. 

Quantum discord has been studied extensively over the 
past few years, from operational definitions 3J, to wit- 
nesses IIHS] and analytical expressions for specific sets 
of states [THT^. However, it is generally difficult to cal- 
culate as it requires both full state knowledge and an 
optimization over all projective measurements of a sub- 
system. Even in the simple case of two qubits, a general 
closed form expression does not exist. The measure of 
Geometric Quantum Discord (GQD) was introduced [5] 
as a simple geometric measure of the distance from a 
given state to the closest classical (zero discord) state. 

DQCl is a model of mixed state quantum computa- 
tion [13] that contains limited entanglement [2] yet is 
thought to outperform classical methods. It does gener- 
ate non-classical correlations as measured by the quan- 
tum discord [6l[Tll[l6], which, in hind sight, is not sur- 
prising since almost all quantum states have non-zero dis- 
cord [yi\. Quantum discord has indeed been witnessed 
in the final states of DQCl computations ^6], even at 
very small values in highly mixed quantum states [B] . To 
further explore the role of quantum correlations in the 
DQCl model and beyond, their quantification becomes 
a necessary pursuit. 

In this work, we derive an analytical expression for 
the geometric discord in the final state of a DQCl com- 
putation of arbitrary dimension, in terms of quantities 
that can be efficiently estimated on a DQCl computer. 



We then demonstrate its experimental evaluation in a 
four-qubit implementation of a DQCl algorithm. Addi- 
tionally, we derive a simple expression for the quantum 
discord in the special case of a two-qubit DQCl-state. 

I. INTRODUCTION 

Quantum discord is the most well-known measure of 
non-classical correlations and is thought to differentiate 
quantum and classical systems [11 12] . It is defined as the 
minimum difference between two classically equivalent 
formulations of the mutual information between subsys- 
tems A and B, I[A : B) := H{A) + H{B) - H{A, B) and 
J{A:B) := H{B) - H{B\A), where H{x) is the Shan- 
non entropy when x is described by a classical prob- 
ability distribution, and is the von Neumann entropy 
when X describes a quantum system. In the quantum 
case the conditional entropy depends on the measure- 
ment basis, H{Uk}{B\A) = Y.kPkH{pB\k), where {11^} is 
a complete set of orthonormal projectors on "H^ such that 
life = I, pfe is the probability of observing outcome k 
on system A, and pB\k ■= TrA[(nfe(g)lB)p(nfc(g)lB)]/pfe is 
the state of system B conditional on the measurement of 
system A returning measurement outcome k. The Quan- 
tum Discord, D{A:B), thus reduces to 

D{A: B) = H{pA)-H{p)+min{n,}Y.PkH^PB\k). (1) 

k 

where pA — TrB(p) is the reduced density matrix of sys- 
tem A. In order to calculate this quantity, full state 
knowledge is required in addition to a minimization over 
all possible projective measurements on subsystem A. 

While there is no analytical expression for the quan- 
tum discord of a general state, the discord for various 
sets of two-qubit states have been found [3 [9l [HI [12] . 
In addition, the quantum discord for all two-qubit states 
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FIG. 1. The circuit diagram above is for the DQCl model 
of computation where the first qubit is in mixed state with 
polarization a along az- At the conclusion of the circuit, 
measurements of (ctj;) and (cry) result in the value of the real 
and imaginary components of the trace of the unitary, scaled 
by the value of the polarization of the first qubit. 

has been reduced to the problem of solving a set of tran- 
scendental equations [TU]. One of the results we report 
here is an expression for the QD in the final state of a 
two-qubit DQCl algorithm that can be calculated with 
only the outcome of the DQCl algorithm (trace of the 
unitary) and the initial polarization of the top register. 

Motivated by the difficultly in computing the QD, Da- 
kic, Vedral and Brukner ^ proposed the GQD, which is 
defined as 

i^^(p)= min ||p-xf , (2) 

where is the set of all zero discord states {D{A : B) = 
0.) These can be written as x = J2jPj\j){j\ ^ pf ■ 
The quantity \\p — xiP = Tr(p — x)^ is the square of 
the Hilbert-Schmidt norm of Hermitian operators. Note 
that the GQD, like the original quantum discord, is not 
symmetric in the subsystems, and Dq{p) is defined as 
minxeof Hp " x\\^ where fl^ = 'EjPjPf ® IjXjI- 

Shortly after the proposal of this measure of quantum 
correlations, Luo and Fu [TH] showed that Equation Q 
is equivalent to the minimization 

D^ip)^rnm\\p-Il^{p)r, (3) 

where = {11^} is a projective measurement on sys- 
tem A, and U^{p) = Eki^t ® ^B)p(n^ (g) Ib)- The 
minimization can be performed analytically for arbitrary 
(2 X 2)-dimensional [5] and (2 x (i)-dimensional [TS] states. 
The resulting expressions are in terms of quantities that 
are not efficiently experimentally accessible. However, a 
method for measuring a tight lower bound on the GQD 
has been recently proposed [20]) requiring only a con- 
stant number of measurements on up to four copies of 
the state. In this report we provide an analytical expres- 
sion for the GQD of (2 x o?) -dimensional DQCl-states, 
and describe a DQCl-algorithm to efficiently estimate it. 



The DQCl model of computation is an example of 
a mixed state quantum information processor |13j . As 
shown in FigurejT] this model has access to only one qubit 
with non-zero polarization, accompanied by a register of 
n maximally mixed qubits. The computation consists of 
two gates: a Hadamard on the top, single qubit register, 
followed by an n-qubit unitary on the bottom register 
that is controlled by the state of the top qubit. The final 
state of this simple algorithm is what we call the DQCl- 
state and can be written as 

1 /^/»" aUl\ 

Expectation value measurements of the Pauli matrices 
and (Jy yield the real and imaginary values of the trace 
of the unitary: ^Re(Tr(C/„)) and ^Im(Tr(C/„)). Since 
calculating the trace of a unitary does not have a known 
efficient classical algorithm, a DQCl computer is believed 
to be more powerful than its classical counterpart. 

Despite its apparent improvement over classical meth- 
ods, DQCl is known to have zero bipartite entanglement 
between the top and bottom registers and can only have a 
small amount of entanglement across any other bipartite 
splitting |14| . It can, however, have quantum correla- 
tions characterized by the quantum discord [15]. This 
has been confirmed in two experiments to date: mea- 
surement of the discord in a two-qubit optics setup using 
full state tomography ^16j , and a four-qubit NMR imple- 
mentation that witnesses discord with a small number of 
experiments 6 . 

II. ANALYTICAL EXPRESSION FOR GQD 

In order to find an expression for the geometric discord 
in a DQCl-state, let us write Eqn. ^ as 

Daip) = min (Tr(p2)-2Tr(pn^(p)) + Tr(n-4(p)2)). (4) 

The first term, the purity of the state of the total system, 
is invariant under unitary transformations. Therefore it 
can be calculated for the initial state, and only depends 
on the initial polarization of the top register: 

Tr(pl3Qci) = ^^^TT- (5) 

In order to calculate the other two terms, consider 
parameterizing the measurement on the top qubit as 
n± = |V'±)(^±I, where \^+) = a\Q) + be"t'\l) and = 
6|0) - ae'-^ll); a and 6 = Vl - are real and e [0, 1]; 
and € [— ^, |]. Thus, the state after measurement is 



PDQCl ® /„) -I- (8) /„) PDQCl (IV'-XV'-I ® In) 

~ 2"+i 2aa^b^iU + e^">'W) /„ - aab{a^ - b^){e-"''U + e^'>'W) ) ' 
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For the DQCl-state, the terms Tt{pII'^{pdqci)) and 
Tr(n^(p£iQCi)^) are equivalent, and evaluate to 
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cos(2(0~0,)), 



where {0j}j=i...2" are the eigenphases of U (i.e. U = 
E^Iie*^'l^j)(^il)> and we have used iTr(e2^'^C/t2 + 
e-2i<^[/2) = 5^^. cos(2((/) - Oj)). Thus, we are left with 
the task of minimizing 

ff(a,0;{^j}) = \\PDQCI~^^{PDQCIW 
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l+2a2„2(^_„2) 
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cos(2(0-0,)), 



over the measurement parameters a £ [0, 1] and € 
[~f ' f ]■ Examining first and second order derivatives 
of 17(0, (j)) with respect to a and 0, we find the optimal 
measurement parameters to be 
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and 



- /E,sin(2(?,) 

Po = arctan =^^^ — — 

2 VE,cos(26lj) 
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arg(Tr((72)), 



which reduces the GQD to 



^g(pdqci)- (|)"^[1-T2] 



(7) 
(8) 
(9) 

(10) 



where T2 = |Tr(C/2)|/2", and can be evaluated with a 
DQCl-algorithm using back-to-back applications of the 
control- J7. The depth of the circuit that evaluates T2 is 
at most double the one that evaluates Tr([/), and there- 
fore has the same efficiency of evaluation as the standard 
DQCl algorithm. Worth noting here that aq = l/\/2 
indicates that the optimal measurement is always in the 
transverse (X-Y) plane of the Bloch sphere, independent 
of the implemented unitary. Also note that Eqn. ([9| al- 
lows for the ability to experimentally evaluate the op- 
timal measurement parameter (jjo using the same data 
collected for the evaluation of T2. 

The GQD does not depend on the eigenvectors of the 
unitary, but rather the distribution of eigenphases. This 
gives rise to classes of unitaries that generate the same 
GQD in a DQCl circuit. For instance, reproducing a 
result from Ref. [S] , the set of unitaries that generate zero 
GQD (and hence, zero QD) will have T2 = 1, implying 
Tr([/2) = e'^2" for some This is possible if and only 
if the eigenphases of U are {dj; 26 j ^ ± 27r V j}. That 
is to say U = e^^^'^A^ where A is a binary observable 

(A2=J). _ 

For a unitary drawn randomly according to the Haar 
measure, its eigenphases are randomly distributed over 



the unit circle |21| . Thus, for large n, T2 approaches 
zero and Dq{pdqci) approaches its maximum value of 
(^)2 2~", i.e. the GQD of a DQCl-state decreases expo- 
nentially with the number of qubits. This is in contrast 
to the average QD of a DQCl-state, which was shown in 
Ref. [IS] to be independent of the number of qubits for 
large n. 



III. EXPERIMENTAL MEASUREMENT OF 
THE GEOMETRIC DISCORD 

We measure the quantum correlations, as quantified 
by the GQD, in a liquid-state nuclear magnetic res- 
onance implementation of the DQCl algorithm. The 
qubits in NMR are ensembles of spin-1/2 nuclei. In 
this experiment, we use the four carbon-13 nuclei in 
the molecule trans-crotonic acid (molecular information 
and Hamiltonian parameters can be found in Ref. |22).) 
The experiment is implemented in a 16.7 T magnetic 
field Bruker Avance spectrometer, where, for carbon-13 
nuclei (7 = 6.728284 x 10^ rad T"^s-i) at room tem- 
perature, the initial polarization of the thermal state is 
a — 1.4x 10~^. The spins are manipulated with radio fre- 
quency pulses numerically generated using the GRAPE 
algorithm [^Sl [H]. The pulses are designed to have a 
fidelity of no less than 0.998 and are adjusted for non- 
linearities in the pulse generation and transmission by 
placing a pickup coil at the location of the sample and 
running a feedback loop to iteratively adjust the pulse 
shape for optimal transmission. 

For the particular instance of the DQCl circuit in this 
experiment, we choose to implement one of the unitary 
matrices used for approximating the Jones polynomial 
for braids with four strands [2?. This is a problem that 
completely encapsulates the power of the DQCl model 
and is of interest in many fields of physics and math. 
The unitary has the form U = diag(c, c, c?, 1, c, d, 1, 1), 
where c = -(e-*3''/5)4 and d = {e''^''/^f. The DQCl- 
state, after the application of the control-?/, is known 
to have nonzero discord as witnessed in a recent exper- 
iment [6]. Using Eqn. (10 1, we calculate the expected 
value of the GQD to be 0.0266 a^. As shown in Fig- 



ure [2j this is achieved for an optimal measurement with 
flo — l/\/2 and (/)o — 0.116 rad. Also shown in Fig- 
ure [2jii) is the experimental value of GQD in the DQCl- 
state, which is found to be (0.0260 ± 0.0004) by mea- 
suring the outcome of a DQCl circuit with back-to-back 
applications of control-?/. For our experimental polar- 
ization [a — 1.4 X 10^^), the GQD in the final state is 
(5.10 ±0.08) X 10-12. 



IV. THE QUANTUM DISCORD OF 
TWO-QUBIT DQCl-STATES 

In the case where the bottom register is a single qubit, 
we report an analytical expression for the QD of the 
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FIG. 2. Shown are (i) the contour plot of the geometric dis- 
tance, p(a, = ||pDQci-n'^(pi5QCi)|P, normalized by Q^, as 
a function of the measurement parameters a and <j), for poQCi 
the output of the four-qubit DQCl algorithm described in the 
text, and (ii) the geometric distance at the optimal measure- 
ment axis a = ao = l/\/2. The dashed lines indicate the pa- 
rameters that correspond to the optimal measurement. The 
experimental data point for the geometric discord of poQCi 
is shown on plot (ii) with error bars propagated from experi- 
mental uncertainties and the spectral fit. 



DQCl-state, thereby contributing to the sets of two- 
qubit states for which the QD has been analyticaUy 
solved. Commencing from Eqn. ([T]), a minimization 
must be performed over the conditional entropy term, 
J2kPkH{pB\k), which, in the case of a two-qubit DQCl- 
state, reduces to 



3 = i 
2 



ni/A + x,)il/4-x,) 



where Xj — -^o>-a\J\ — cos((/) — p± = 1/2 ± Xj, 
and a, 6, and parameterize the measurement on the top 
qubit (see the GQD case above). The minimization with 



respect to the measurement parameters a and results 
in an optimal measurement characterized by 

ao--^ and ^,= - + ^^. 

This allows us to write the analytical expression for the 
QD of a two-qubit DQCl-state as 



D{A : B) =H2 



1 — a Ti 



-Ho 



1 - a 



l-T^log 




(11) 



where i?2(') is the binary entropy, and Ti — |Tr(C/)|/4, 
which can be measured directly with a DQCl-algorithm. 

Hence, the quantum discord of a two-qubit DQCl- 
state is zero for |Tr(?7)p = or 4, and maximum for 
|Tr(C/)|2 ^ 2. RecaU, from Eqn. that the GQD is 

zero for \Tr{U^)\ = 1/2" and maximum for |Tr(t72)| = 0. 
Consider a general single-qubit unitary, expressed as 
U = e'^^Rf{'&), where Rf{^) ■= exp[— it9f • (7/2] is a ro- 
tation by "i? about the unit vector in three dimensions 
f = (rx,ry,rz)^ o — {<yx,<^y,'^z) are the Pauli matri- 
ces, and ^ a global phase. It follows from (f • a)^ = I 
that Tt{U) = 2e*« cos(i?/2) and Tt{U^) = 2e*2« cos(t?). 
Thus, the QD and GQD are simultaneously maximum 
for I? = {2k + l)7r/2, and simultaneously zero for ■& = kw, 
independent of the rotation axis f. Examples of uni- 
taries that can be written as tt rotations include bit-flip 
(or NOT) and Hadamard. 



V. CONCLUSION 

In this article we provided analytical expressions for 
the quantum correlations present in a DQCl-state in 
terms of experimentally accessible quantities. The ge- 
ometric quantum discord can be computed for (2 x d)- 
dimensional DQCl-states by implementing a DQCl al- 
gorithm with a second application of the controlled uni- 
tary. The algorithm was experimentally demonstrated 
for a four-qubit liquid-state NMR implementation of the 
DQCl model. We also showed, that for the special case 
of a two-qubit system, the quantum discord can be cal- 
culated using outcome of the DQCl-circuit (the trace of 
the unitary). 

For a typical DQCl-state, where the unitary is chosen 
uniformly at random according to the Haar measure, we 
found that the geometric measure of discord scales with 
2"", in contrast to the quantum discord, which is known 
to be independent of n (for large n) . This result further 
supports the suspicion that the apparent speedup exhib- 
ited by (the dynamics of) DQCl is not captured by the 
geometric measure of quantum discord at the conclusion 
of the algorithm. 
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